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Abstract. We apply the idea that dense stellar systems in the central region of galaxies are 
formed via globular cluster mergers to the formation of the nuclear star cluster of the Milky 
Way, where a massive black hole (Sgr A*) is present. Our high precision N- body simula- 
tions show a good fit to the observational characteristics of the Milky Way nuclear cluster, 
giving further reliability to the so called migratory model for the formation of compact 
systems in the inner galaxy regions 
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1. Introduction 

Compact Massive Objects (CMOs) are almost 
ubiquitously present in the galactic centers, in 
form of massive or supermassive black holes 
(MBHs, SMBHs) or in a more dilute form, i.e. 
resolved stellar nuclei or Nuclear Star Clusters 
(NSCs). There is a direct correlation between 
compactness of the CMO and the host galaxy 
luminosity: the brighter the host the denser the 
CMO, till the limit o f BH i nfinite density (see, 
e.g JVolonteri et al. I ( 120081) "). 

An interpretation of the modes of for- 
mation of CMOs in galaxies is surely a 
modern topic, which have not yet received 
an adequate self-consistent explanation frame. 
Preliminarily, we must say that a common ori- 
gin for the CMOs in galaxies, with a sec- 
ond order difference in their compactness and 
mass as due to properties of the host galaxy, 
is an appealing view, but faces with the ev- 
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idence, for example, that some galaxies con - 
tain both an NSC and a SMBH ( iBoker 1120081) . 
A well-known example is the Milky Way 
(MW), where a 4 x 10 6 M black hole co- 
exists with an NSC ~ 4 times more massive 
(Mjvsc ~ 1-5 x 10 7 Mq); the same coexis- 
tence has been found recently in various other 
cases. Anyway, in spite of these peculiarities, 
it is surely intriguing the p ossibility, raised 
bv lCapuzzo-Dolcetta I d 19931) and furtherly de- 
veloped by him and collaborators, that CMOs 
have formed by stellar aggregations like mas- 
sive globular clusters (GCs) brought to the host 
galaxy center via dynamical braking (dynami- 
cal friction). The infalling and orbitally decay- 
ing GCs eventually merge in the innermost re- 
gion of the host galaxy where they form an 
actual Super Star Cluster in times and modes 
that are, partly, governed by the galaxy struc- 
ture, i.e. geometrical shape, steepness of mass 
density distribution as well as total mass. If the 
two phenomena of dynamical friction braking, 
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Fig. 1. Projected density profile at the end of the 
merger simulation. The solid line gives the best- 
fitting model to the entire system (galaxy+NSC), 
and the dashed curve gives the fit to the density pro- 
file of the galaxy (from Aetall2). 
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Fig. 2. Post-infall evolution of the ,/V-body model. 
The plot shows the break (core) radius as a function 
of time (in units of the relaxation time at the BH 
influence radius). The solid line is the best-fit expo- 
nential (from Aetall2). 



yielding to a central mass accumulation, and 
tidal disruption, working in the opposite di- 
rection, are the leading mechanisms, the host 
galaxy luminosity (mass) correlation with the 
CMO compactness has a qualitative explana- 
tion. Actually, fainter galaxies show a higher 
central density in phase space, thus enhancing 
the effect of dynamical friction on GCs and 
favoring their accumulation to grow a NSC, 
while bright and very bright galaxies borrow, 
usually, MBHs and SMBHs in their center 
which may act as destroyers of incoming GCs 
inibiting the local NSC formation. 

2. The formation of the MW Nuclear 
Star Cluster 

The possibility of cumulation of a signifi- 
cant amount of mass in the innermost re- 
gions of galaxies i n form of infalli n g GC s 
is an old idea by iTremaineT et al. I (I1975I) . 
For some time, the dynamical friction effect 
was considered important just for very mas- 
sive GCs, u ntil its role, enh anced in triax- 
ial galaxiesj Pesce et aiTl 1992b and in cuspy 



il galaxies (fesceetal. iwz) ana in cuspy 
alaxies ( Vicari & Capuzzo-Dolcetta 2005; 



lArca-S edda & Capuzzo-Dolcetta 1 1201 2l) . was 
clearly shown to be relevant for a wide range of 
G C structural parameters and initial conditions 
bv ICapuzzo-Dolcetta 1 (11993b and later papers. 
This framework of GC infalls provides possi- 
ble answers on both the origin of local black 
hole mass feeding and on the modes of for- 
mation of the observed dense stellar systems 
in both dwarf ellipticals and spirals. Actually, 
these systems have characteristics similar to 
those of GCs but on a magnified scale. This mi- 
gratory explanation is alternative to that claim- 
ing a funneling of molecular gas toward the 
inner part of galaxies and subs equent local 
star fo rmation (see for instance lLoose et al. I 
(Il982l) '). The two hypotheses have been dif- 
ferently tested by various authors; the dissi- 
pationless hypothesis (GC migratory model) 
has also been quantitatively well tested by spe- 
cific A^-body simulations while the dissipative 
(gaseous) origin has not yet been thoroughly 
investigated. 

It seems, so, logical performing a test of 
the migratory model for the origin of galactic 
central dense stellar systems in the best known 
environment: the Milky Way (MW), where, 
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around the ~ 4 x 10 6 M central black hole 
there is a NSC of about 10 7 M Q whose pho- 
tometric structure as well as kinematical char- 
acteristics are studied at a level unreachable in 
external galaxies. As a matter of fact, the MW 
NSC is, so far, the only one resolved in its stel- 
lar content in spite of its hig h density. As a 
conseq uence, in a recent paper lAntonini et akl 
d2012l) (hereafter Ae tall 2) tackled the problem 
of analyzing whether the migratory hypothesis 
for the formation is working also in the specific 
context of the MW, where, indeed, the dynam- 
ical role of the black hole is likely to be signif- 
icant. 

The t observational constraints are: 
the mass of the NSC (~ 1.5 x 10 7 M ); ii) 
its mass density profile, p(r), showing a core in 
the distribution of the late type stars of 0.5 pc 
size, and decreasing further out as p oc r~ 18 up 
to r ^ 30 pc. Out of 30 pc there is a large nu- 
clear stellar and molecular disk of about same 
radius; in) some kinematic data; iv) the NSC 
luminosity function. 

We showed how the orbital decay of 12 
tidally truncated massive (~ 10 6 M ) GCs ac- 
tually leads, in a sufficiently short time, to the 
growth and stabilization of a dense stellar clus- 
ter around the Sgr A* BH. We followed the in- 
fall and merger of the GCs in two sets of ini- 
tial conditions (simultaneous and consecutive 
infalls) and were able to study the density pro- 
file shape as well as the distribution in veloc- 
ity space. Moreover, we studied also the post 
merger evolution, extending the A^-body sim- 
ulations after the final merger for a time that 
corresponds to about 10 Gyr, after scaling to 
the MW. The absence of a clear formation of a 
Bahcall-Wolf cusp in the NSC is not surprising 
because the expected time scale for its growth 
is the 2-body relaxation time at the black hole 
influence radiu s is in the range 20 - 30 Gyr 
(lMerrittll2010h . 

3. Results 

Taking into account the constraints imposed by 
the observations of the MW NSC, our results 
are conveniently summarized as follows: 

- The stellar system resulting from the con- 
secutive mergers has a density that falls off 



as r 2 and a core of radius ~ 1 pc. These 
properties are similar to those observed in 
the MW NSC (see Fig. Q). 

- The morphology of the NSC evolved dur- 
ing the infalls, from a strongly triaxial 
shape toward a more oblate/axisymmetric 
shape. Kinematically, the final system 
is characterized by a mild tangential 
anisotropy within the inner 30 pc and a low 
degree of rotation. 

- The effect of gravitational encounters on 
the evolution of the NSC was investigated 
continuing A^-body integrations after the 
last inspiral event. The NSC core shrinks 
by roughly a factor of two in 10 Gyr as 
the stellar density starts evolving toward 
a Bahcall-Wolf cusp (see Fig. [2]). This fi- 
nal core size is essentially identical to the 
size of the core observed at the center of 
the Milky Way. The density profile outside 
the core remained nearly unchanged dur- 
ing this evolution. Gravitational encounters 
also make the NSC to evolve toward spher- 
ical symmetry in configuration and veloc- 
ity space. 

- Since GCs are ancient objects with ages 
1013 Gyr, in the merger model a large frac- 
tion of the NSC mass is expected to be in 
old stars. Using stellar population models, 
we showed that the observed LF of the MW 
NSC is consistent with a star formation his- 
tory in which a large fraction (about 1/2) of 
the mass consists of old (10 Gyr) stars and 
the remainder is from continuous star for- 
mation (see Fig . [3j . 

4. Conclusions 

The migratory-merger model for the formation 
of dense stellar systems in the centers of ellip- 
tical galaxies and of superdense Nuclear Star 
Clusters in spirals has been applied to the case 
of the Milky Way Nuclear Star Cluster, where 
a massive black hole dominates the innermost 
potential. Our A^-body simulations show that 
a NSC with the same spatial and kinemati- 
cal characteristics of the MW NSC is formed 
by inspiralling of a limited number (~ 10 — 
15) massive globular clusters in a time short 
respect to the age of the galaxy. Moreover, 
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Fig. 3. Left panels display star formation scenarios in which, from top to bottom: the entire mass is built 
up by ancient stars, brought into the Galactic center by infalling GCs; half of the mass is contributed by 
old stars and half due to continuous star formation; at the mass contributed by continuous star formation 
and GCs we add stars formed during a starburst episode occurred at ~ 1 Gyr; all the mass is contributed by 
continuous star formation. Right panels compare the observed LF of the Galactic center (light lines) with 
the LFs (heavy lines) resulting from the different star formation scenarios assuming solar metallicity and 
canonical mass-loss rates. Note that the data are much more than 50% incomplete for the faintest few bins. 
The models have not been scaled for mass, but rather have been rescaled along the vertical axis to match 
the number counts in the K = 1 1.0 bin (from Aetall2). 



the observed luminosity function of the MW 
NSC is compatible with a younger population 
emebedded in that coming from the merger ob- 
ject. 
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